Research Article
Enrichment Analysis of Differentially Expressed Genes during Endosperm Grouting Periods in Non-waxy and Waxy Foxtail Millets ![](https://genbreedpublisher.com/templates/images/OA.gif)
2 Crop institute, Inner Mongolia Academy of Agriculture and Animal Husbandry, Hohhot, 010031, China
* These authors contributed equally to this work
![](https://genbreedpublisher.com/templates/images/email1.bmp)
![](https://genbreedpublisher.com/templates/images/email2.bmp)
Molecular Plant Breeding, 2020, Vol. 11, No. 12 doi: 10.5376/mpb.2020.11.0012
Received: 18 Jun., 2020 Accepted: 28 Aug., 2020 Published: 30 Sep., 2020
Li Q., Li T., Guo S.H., Bai Y.T., and Li X.C., 2020, Enrichment analysis of differentially expressed genes during endosperm grouting periods in non-waxy and waxy foxtail millets, Molecular Plant Breeding, 11(12): 1-8 (doi: 10.5376/mpb.2020.11.0012)
The content of amylose and amylopectin is an important trait affecting cooking and eating quality in millet. During process of kernel grouting, the synthesis of starch involves different pathways and components in non-waxy and waxy millets. Immature grouting grains of waxy millet ‘Gonggu 68’ and non-waxy millet ‘Chigu 4’ (grouting period 1 and 5 days) were used to analyze their transcriptome sequences by using Illumina Hiseq4000. The results showed that: (1) GBSSⅠenzymes activity of waxy and non-waxy millet was low-high-low. There were some differences between the two activities. 665 up-regulated differentially expressed genes were screened on day 5 and day 1 during grouting period in waxy cultivar ‘Gonggu 68’, there were 431 more up-regulated genes than down-regulated genes. There were 97 more up-regulated genes than down-regulated genes in non-waxy cultivar ‘Chigu 4’ on day 5 and day 1 in grouting period. (2) In the A2-VS-A1 waxy comparison group, the differential genes were mainly GO enriched in 7 functions such as the seed oil body biogenic function, 17-β-ketosteroid reductase activity function and so on. it was mainly enriched in biological processes and molecular functions. In the B2-VS-B1 non-waxy comparison group, the differential genes were mainly GO enriched in 8 functions such as the light harvesting in photosystem I function and the pigment binding function for non-waxy millet. (3) Differentially expressed genes were mainly KEGG enriched in caffeine metabolism pathway, linoleic acid metabolism pathway, anthocyanin biosynthesis pathway, aflatoxin biosynthesis pathway in waxy A2-VS-A1, but which were mainly KEGG enriched in the synergy-antenna protein pathway, the linoleic acid metabolic pathway, the caffeine metabolic pathway, the brassinosteroid biopathway in non-waxy B2-VS-B1. These two comparative groups were enriched Caffeine metabolism pathway and linoleic acid metabolism pathway appeared in the process. (4) Three (SSII-3, PHO1, AS) and four (PHO1-1, AS, AGP16, WAXY) genes with significant differences and related to waxy and non-waxy millet were screened. With Actin (Si001873) as the internal reference gene, the above seven differentially expressed genes were verified by qRT-PCR, which was consistent with the transcriptome results, indicating that the differentially expressed genes were related to waxy or non-waxy endosperm.
Setaria italica (L.) Beauv. belongs to the Poaceae Barnhart, also known as Daliang or Xiaomi in China. Setaria italica (L.) Beauv. is the earliest cultivated grain in primitive agriculture in northern China. For thousands of years, Setaria italica (L.) Beauv. has been widely planted because of its favorable characteristics such as drought tolerance, fast growth and strong adaptability, which is distributed from Xinjiang in the west to the whole Yellow River Basin (You, 1993, Agricultural History of China, (3): 1-13). Setaria italica (L.) Beauv. is rich in protein, fat, carbohydrates, vitamins and minerals and provides six basic nutrients for the human body (Qiao and Wang, 2015, Agricultural Science and Technology and Equipment, (11): 41-42). Due to different planting areas and varieties, the content of each nutrient may be different (Xue et al., 2008). In recent years, people pay more and more attention to eating habits and dietary structure, and Setaria italica (L.) Beauv. is favored by more and more people. Therefore, it has gradually become an indispensable part of the diet and been one of the most popular food crops in China.
In higher plants, starch can be divided into amylose and amylopectin according to different functions and structures. According to amylose content, it can be divided into non-waxy millet and waxy millet. The amylose content of non-waxy millet is >17.0%, and that of waxy millet is 0~3.5% (Nakayama et al., 1998). Waxy gene encodes granule-bound starch synthase (GBSSⅠ), which is a key enzyme for amylose synthesis (Yang et al., 2007) and widely exists in rice, corn and other food crops (Gao et al., 2006). Deletion or mutation of waxy gene will affect the activity of GBSSⅠ, which in turn affects the synthesis of amylose (Zhang et al., 2009). The content of amylose and amylopectin is an important trait that affects the cooking and eating quality of millet (Jia et al., 2010). Waxy millet is deeply loved by consumers because of its rich nutrition and diverse tastes. Therefore, it is necessary to study the mechanism of waxy millet and explore genes related to waxy phenotype. At present, there are no reports about the use of transcription sequencing technology to discover the differentially expressed genes during grouting periods in millets in China.
In this study, ‘Chigu 4’ (non-waxy millet) and ‘Gonggu 68’ (waxy millet) were selected to determine the activity of granule-bound starch synthase (GBSSⅠ) during grouting periods. The grain samples corresponding to the time points with the highest and lowest enzyme activity were selected for transcriptome sequencing, screening, verification and enrichment analysis of differentially expressed genes.
1 Results and Analysis
1.1 The determination of GBSSⅠ enzyme activity
The enzyme activities of ‘Chigu 4’ and ‘Gonggu 68’ were lowest on day 1 during grouting periods, which were 200 and 300 IU respectively (Figure 1), and the GBSSⅠ activities of the two varieties were highest on day 5, which were 660 and 560 IU respectively.
![]() Figure 1 Dynamic changes of granule-bound starch synthase (GBSSⅠ) activity |
1.2 The analysis of differentially expressed genes (DEGs)
From the number distribution of differentially expressed genes in waxy and non-waxy millet cultivar (Figure 2), it can be seen that 665 up-regulated and 234 down-regulated differentially expressed genes were screened from the A2-VS-A1 comparison group on day 5 and day 1 during grouting period of waxy cultivar ‘Gonggu 68’, there were 431 more up-regulated genes than down-regulated genes. While 1 586 up-regulated and 1 489 down-regulated differentially expressed genes were screened from the B2-VS-B1 comparison group on day 5 and day 1 during grouting period of non-waxy cultivar ‘Chigu 4’, there were 97 more up-regulated genes than down-regulated genes.
![]() Figure 2 Number distribution of differentially expressed genes in waxy and non-waxy millet cultivar |
1.3 GO enrichment analysis of differentially expressed genes
In the A2-VS-A1 waxy comparison group, differential expressed genes (DEGs) were enriched in many pathways (Figure 3). The enrichment factors that were larger than 0.1 included the following six pathways, namely seed oilbody biogenesis function, 17-beta-ketosteroid reductase activity function, 17-beta-hydroxysteroid dehydrogenase (NADP+) activity function, killing of cells of other organism function, central vacuole function, lipid storage function. Among them, seed oilbody biogenesis function, killing of cells of other organism function and lipid storage function were biological processes, while 17-beta-ketosteroid reductase activity function and 17-beta-hydroxysteroid dehydrogenase (NADP+) activity function were molecular functions. In the B2-VS-B1 non-waxy comparison group, the enrichment factors that were larger than 0.15 also included the following six pathways, namely light harvesting in photosystem I function, pigment binding function, chlorophyll binding function, pectinesterase activity function, pectinesterase inhibitor activity function, plastoglobule function, cellular transition metal ion homeostasis function, transition metal ion binding function and so on. Among them, light harvesting in photosystem I and cellular transition metal ion homeostasis function were biological processes, while pigment binding, chlorophyll binding, pectinesterase activity, pectinesterase inhibitor activity, plastoglobule function, and transition metal ion binding and extracellular regional function were molecular functions.
![]() Figure 3 GO enrichment analysis of differentially expressed genes |
1.4 KEGG enrichment analysis of differentially expressed genes
In the A2-VS-A1 waxy comparison group (Figure 4), the enrichment factors of the differentially expressed genes that were larger than 0.06 included caffeine metabolism pathway, linoleic acid metabolism pathway, anthocyanin biosynthesis pathway, aflatoxin biosynthesis pathway and so on. In the B2-VS-B1 non-waxy comparison group, the enrichment factors of the differentially expressed genes that were larger than 0.06 included photosynthesis-antenna protein pathway, linoleic acid metabolism pathway, caffeine metabolism pathway and brassinosteroid biosynthesis pathway.
![]() Figure 4 KEGG enrichment analysis of differentially expressed genes |
1.5 Screening and verification of differentially expressed genes
Using |log2FC|>1, p<0.05 as the screening conditions, the differential expressed genes in the two comparison groups were screened. There were 899 and 6 150 differential expressed genes in A2-VS-A1 and B2-VS-B1, respectively. According to the gene function annotation, 36 genes related to the waxy and non-waxy millet were screened, among which, 12 was in A2-VS-A1 group and 24 was in B2-VS-B1 group. The results (Table 1) showed that 3 differential genes with FDR<0.01 and high expression (TPM), related to the waxyand non-waxy millet, were screened in A2-VS-A1 comparison group, including SSII-3, PHO1 and AS, while in B2-VS-B1 group, 4 differential genes were screened, including PHO1-1, AS, AGP16 and WAXY.
![]() Table 1 Differential genes and their expression quantity |
According to the screening results of differential genes, the differential genes were verified by real-time fluorescence quantitative qRT-PCR with Actin (Si001873) as the internal reference gene. The total RNA of the sample was extracted, and the quality and purity of the sample were detected by 1.5% agarose gel electrophoresis (Figure 5). The differential genes were performed t test at the 0.05 level by SPSS 19.0 (Figure 6).
![]() Figure 5 Agarose gel electrophoresis of total RNA extraction |
![]() Figure 6 New ICT based fertility management model in private dairy farm India as well as abroad |
In A2-VS-A1 waxy comparison group (Figure 6), A1 was the control group, and the expression of differential gene SSⅡ-3 and PHO1 inA2 was lower than that in A1, showing down-regulation expression; The differential gene AS showed up-regulation expression. In B2-VS-B1 non-waxy comparison group, B1 was the control group, and the expression of differential genes PHO1-1, AGP16 and WAXY inB2 were lower than that in B1, showing down-regulation expression, but the expression of differential gene AS in B2 was higher than that in B1, showing up-regulation expression.
2 Discussion
The differentially expressed genes can be identified by comparing the high-throughput transcription and sequencing data of different tissues. New transcripts, variable splicing and gene structure variation can be identified and discovered with the help of a large number of gene annotation data (Qi et al., 2011). With the application of bioinformatics in the field of life science, transcriptome sequencing has become an important research method of transcriptome. Transcriptome sequencing technology is relatively mature, low-cost and widely used, which comprehensively and quickly reflects the transcriptome data of a certain tissue or cell at a specific time and specific location (Zhou et al., 2019). The researchers obtained the anthocyanin metabolic structure gene and regulatory gene sequence of Hyacinthus orientalis L. variety ‘Black Baccara’ through transcriptome sequencing. The expression patterns of ‘Black Baccara’ in white, red, blue and black were analyzed by qRT-PCR. It was found that HoF3H, HoMYB4, HoMYB6 and HoMYB7 were four key genes affecting anthocyanin accumulation of Hyacinthus orientalis L. (Li et al., 2020, http://kns.cnki.net/kcms/detail/46.1068.S.20200226.1514.004.html). Wang et al. (2018) used high-throughput RNA-seq technology to sequence the transcriptome of drought-stressed millet, and the samples treated 6 d were compared with that treated 6 hours. It was found that 206 genes were up-regulated and 611 genes were down-regulated. Zu et al. (2015) obtained a WRKY transcription factor SiWRKY36 by sequencing the transcriptome of drought-stressed millet.
In this study, the high-throughput sequencing technology Illumina Hiseq4000 was used to sequence the transcriptome of immature grains at a specific time during the grouting period of non-waxy millet (‘Chigu 4’) and waxy millet (‘Gonggu 68’). There were 431 more up-regulated genes than down-regulated genes in the waxy comparison group A2-VS-A1, while there were 97 more up-regulated genes than down-regulated genes in the non-waxy comparison group B2-VS-B1. The results showed that the number of up-regulated genes was significantly larger than that of down-regulated genes in the early grouting period of waxy millet, and there was less difference between up-regulated genes and down-regulated genes in the early grouting period of non-waxy millet.
In the A2-VS-A1 waxy comparison group, the differential genes were mainly GO enriched in 7 functions such as the seed oil body biogenic function, 17-β-ketosteroid reductase activity function and so on. In the B2-VS-B1 non-waxy comparison group, the differential genes were mainly GO enriched in 8 functions such as the light harvesting in photosystem I function and the pigment binding function for non-waxy millet. The two comparison groups were mainly enriched in biological processes and molecular functions.
In the A2-VS-A1 waxy comparison group, the differential genes were mainly KEGG enriched in caffeine metabolism pathway, linoleic acid metabolism pathway, anthocyanin biosynthesis pathway, aflatoxin biosynthesis pathway and so on. In the B2-VS-B1 non-waxy comparison group, the differential genes were mainly KEGG enriched in photosynthesis-antenna protein pathway, linoleic acid metabolism pathway, caffeine metabolism pathway and brassinosteroid biosynthesis pathway. Both the non-waxy millet and waxy millet included caffeine metabolism pathway and linoleic acid metabolism pathway.
Soluble starch synthase genes include SSI, SSII, and SSIII, which participate in amylopectin synthesis together with branching enzyme genes and debranching enzyme genes (Ohdan et al., 2005; Pandey et al., 2012). The SSII-3 gene is the main gene that controls the gelatinization temperature of rice (Gao et al., 2003, Science in China (Serices C: Life Sciences), 33(6): 481-487). SSII-3 has two genotypes, and different genotypes have different effects on eating quality. Type I is high gelatinization temperature and low amylose, and type II is low gelatinization temperature and high amylose (Tian et al., 2009). In this study, the expression of SSII-3 gene was down-regulated in A2-VS-A1 group, that is, the expression of SSII-3 gene in the middle grouting period of ‘Gonggu 68’ was lower than that in the early grouting period, indicating that the amylopectin synthesis in waxy millet was more active in the early grouting period, which may be an important period to determine the non-waxy and waxy millet and a period that has a great impact on the eating quality of millet.
AS is one of glycosyltransferase genes. AS can accept a variety of compounds with diverse structures as its substrate, which is a rare multifunctional glycosyltransferase with multi-substrate acceptance ability (Hefner et al., 2002). In this study, the expression of AS on day 5 was higher than that on day 1 in the two comparison groups, indicating that the expression of AS increased continuously in the early grouting period, which was more important in the process of grain formation.
Starch phosphorylase plays a dynamic regulatory role in the carbohydrate metabolism of higher plants, that is, the synthesis and degradation of starch (Schneider et al., 1981). In this study, the PHO1 gene showed a down-regulation trend in the two comparison groups, and waxy millet had an obvious down-regulation trend, indicating that the expression of starch phosphorylase gene was higher in non-waxy millet and waxy millet in the early grouting period, which was more important for starch synthesis and degradation.
In the process of starch synthesis, adp-glucose pyrophosphorylase (AGP) catalyzes glucose-1-phosphate (G-1-P) and adenosine triphosphate (ATP) to synthesize adenosine diphosphate glucose (ADPG) and release pyrophosphate (PPi) (Gao and Huang, 1998). The effect of genotype on AGP activity and starch content is simultaneous, and AGP activity has a direct effect on starch content (Zhao et al., 2013).
In this study, the expression of differential gene AGP16 was down-regulated in B2-VS-B1 group, indicating that the expression of differential gene AGP16 in the middle grouting period was lower than that in the early grouting period, and the early grouting period was an important stage for starch formation. AGP has a variety of isozymes, and AGP16 gene may have similar function to AGP gene, which may be its isozyme gene.
Waxy gene controls the amylose synthesis in endosperm and encodes a key enzyme - granule-bound starch synthase (GBSSⅠ) for plant starch synthesis. The hindrance of amylose synthesis in endosperm and the decrease of content is caused by the deletion or mutation of Waxy gene or the decrease of GBSSⅠ enzyme activity (Liu et al., 2010). In this study, the expression of differential gene WAXY was down-regulated in B2-VS-B1 group, indicating that the expression of differential gene WAXY in the middle grouting period of non-waxy millet was lower than that in early grouting period, and the early grouting period of non-waxy millet was an important stage for amylose synthesis. The formation of waxy endosperm is caused by the deletion and mutation of WAXY gene.
In this study, 7 differentially expressed genes in different periods of non-waxy and waxy millet were obtained by transcriptome sequencing, and the differentially expressed genes were verified by qRT-PCR, which provided an experimental basis for millet quality genetic breeding, germplasm identification and molecular marker assisted breeding.
3 Materials and Methods
3.1 Plant materials
The materials of this study were two spring millet varieties introduced from China: ‘Gonggu 68’ (waxy variety, AC was 0, selected by Jilin Academy of Agricultural Sciences) and ‘Chigu 4’ (non-waxy variety, AC was 20.10%, selected by Chifeng Academy of Agricultural Sciences).
3.2 Sampling and enzyme activity determination
In this study, the two experimental materials were designed according to the random block. The plot area was 25 m2 and repeated for 3 times, which was planted on the Inner Mongolia Agricultural Teaching Farm. The 12th day after millet flowering was selected as the first day of grouting period. The sampling time points were day 1, day 3, day 5, day 10, day 15, day 20, day 25 and day 30. After liquid nitrogen quick freezing, transfer to ultra-low temperature refrigerator to save for testing. The enzyme activity was measured by using the GBSSⅠ kit from Suzhou Comin Biotechnology Co., Ltd. According to the level of enzyme activity, comparison groups were set for the highest and lowest points of enzyme activity of waxy millet and non-waxy millet, named A2-A1, B2-B1, and transcriptome sequencing was performed on the corresponding grain materials to compare the expression of differential genes in the same variety at different grouting periods.
3.3 Analysis of differentially expressed genes (DEGs) and gene function and classification annotation
According to the expression of genes in the comparison groups (TPM) in the sequencing results, DEGs were screened, and the genes were annotated with GO and KEGG pathways.
3.4 Real time fluorescence quantitative PCR analysis
The total RNA was isolated from the corresponding grain materials and its quality and purity were determined by agarose gel electrophoresis. Taking Actin (Si001873) as the internal reference gene, the gene sequence primers were designed by PrimerQuest tool (Table 2). The expression of DEGs was calculated by 2-∆∆Ct method and plotted by GraphPad Prism 8.0.1.
![]() Table 2 Candidate genes and their expression levels in non-waxy and waxy comparable group Note: A: Candidate genes and their expression levels in waxy comparison group; B: Candidate genes and their expression levels in non-waxy comparison group |
Authors’ contributions
GSH, LQ, and LT were the designers of the experiments and executors of the research; LQ and LT completed the data analysis and wrote the first draft of the manuscript; BYT and LXC participated in designing the experiments and analyzing the experimental results; GSH was the conceiver and person in charge of the project, directing experimental design, data analysis, manuscript writing and revision. All authors read and approved the final manuscript.
Acknowledgments
This study was jointly funded by the National Natural Science Foundation of China (31560394), the Science and Technology Plan of Inner Mongolia Autonomous Region (201802064) and the Special Fund of Breeding of Animal and Plant Varieties for Transformation of Scientific and Technological Achievements of Inner Mongolia Agricultural University (YZGC2017022).
Gao D.R., Guo W.S., Cheng S.H., Feng C.N., Zhu X.K., and Peng Y.X., 2006, Characteristics of pre-harvest sprouting in wheat with different types of Wx protein deficiency, Zuowu Xuebao (Acta Agronomica Sinica), 32(12): 1917-1919
Gao Z.Y., and Huang D.N., 1998, Starch branching enzyme: The key enzyme of amylopectin biosynthesis in plants, Shengwu Gongcheng Jinzhan (Progress in Biotechnology), 18(3): 26-29
Hefner T., Arend J., Warzecha H., Siems K., and Stöckigt J., 2002, Arbutin synthase, a novel member of the NRD1 glycotransferase family, is a unique multifunctional enzyme converting various natural products and xenobiotics, Bioorg. Med. Chem., 10(6): 1731-1741
https://doi.org/10.1016/S0968-0896(02)00029-9
Jia L.X., Ma N.P., Jing W.W., and Xue L.T., and Zai Y.D., 2010, Determination of amylose and amylopectin in grains of cereals and Legume by dual-wavelength spectrophotometry, Xinjiang Nongye Kexue (Xinjiang Agricultural Sciences), 47(3): 564-568
Liu H.M., Zhao Y., Gu Y., Zhang H.Y., and Huang Y.B., 2010, Characterization of codon usage of waxy genes in several plants, Henongxue Bao (Journal of Nuclear Agricultural Sciences), 24 (3): 476-481
Nakayama H., Afzal M., and Okuno K., 1998, Intraspecific differentiation and geographical distribution of Wx alleles for low amylase content in endosperm of millet, Setaria italic (L.) Beauv., Euphytica, 102(3): 289-293
https://doi.org/10.1023/A:1018344819531
Ohdan T., Francisco P.B. Jr., Sawada T., Hirose T., Terao T., Satoh H., and Nakamura Y., 2005, Expression profiling of genes involved in starch synthesis in sink and source organs of rice, J. Exp. Bot., 56(422): 3229-3244
https://doi.org/10.1093/jxb/eri292
Pandey M.K., Rani N.S., Madhav M.S., Sundaram R.M., Varaprasad G.S., Sivaranjani A.K.P., Bohra A., RamKumar G., and Kumar A., 2012, Different isoforms of starch synthesizing enzymes controlling amylose and amylopectin content in rice (Oryza sativa L.), Biotechnol. Adv., 30(6): 1697-1706
https://doi.org/10.1016/j.biotechadv.2012.08.011
Qi Y.X., Liu Y.B., and Rong W.H., 2011, RNA-Seq and its applications: a new technology for transcriptomics, Yichuan (Hereditas), 33(11): 1191-1202
https://doi.org/10.3724/SP.J.1005.2011.01191
Schneider E.M., Becker J.U., and Volkmann D., 1981, Biochemical properties of potato phosphorylase change with its intracellular localization as revealed by immunoical methods, Planta, 151(2): 124-134
https://doi.org/10.1007/BF00387813
Tian Z., Qian Q., Liu Q., Liu Q.Q., Yan M.X., Liu X.F., Yan C.J., Liu G., Gao Z.Y,. Tang S.Z., Zeng D., Wang Y.H., Yu J.M., Gu M.H., and Li J.Y., 2009, Allelic diversities in rice starch biosynthesis lead to a diverse array of rice eating and cooking qualities, Proc. Natl. Acad. Sci. USA, 106(51): 21760-21765
https://doi.org/10.1073/pnas.0912396106
Wang Q.G., Li Z., Guan Y.A., Liu W., and Pan J.W., 2018, Transcriptome profiling analysis of foxtail millet under drought stress, Shandong Nongye Kexue (Shandong Agricultural Sciences), 50(10): 6-12
Xue Y.Y., Li P., and Lin Q.B., 2008, Research evolution on chemical component and physical character of foxtail millet, Zhongguo Liangyou Xuebao (Journal of the Chinese Cereals and Oils Association), 23(3): 199-203
Yang R.G., Dai Z.M., Wang H.M., Yu J., Li P.H., Li Z.Q., and Wu Z.J., 2007, The study progress in biosynthesis of plant starch, Anhui Nongye Kexue (Anhui Agricultural Science), 35(28): 8791-8793
Zhang M., Li Y.D., Ma J.F., Liu H., Li X.L., Xu J., and Tang L., 2009, Quality analysis on 7 millet cultivars in Heilongjiang area, Shipin Kexue (Food Science), 30(1): 41-43
Zhao F.M., Yue X.W., Yang H.H., Xiao J., Li Y.X., and Wang X.Z., 2013, Isozyme forms and temporal and spatial expression profile of ADP-glucose pyrophosphorylase in developing wheat grain, Xibei Zhiwu Xuebao (Acta Botanica Boreali-Occidentalia Sinica), 33(8): 1528-1534
Zhou W., Hu N., Liu X.T., Wang Y.W., and Suo Y.R., 2019, The RNA-Seq transcriptome analysis of female and male plants leaves in Hippophae rhamnoides L. Subsp. sinensis, Jiyinzuxue Yu Yingyong Shengwuxue (Genomics and Applied Biology), 38(12): 5516-5526
Zu Q.L., Yin L.J., Xu Z.S., Chen M., Zhou Y.B., Li L.C., Ma Y.Z., Min D.H., and Zhang X.H., 2015, Molecular characteristics and functional identification of foxtail millet transcription factor WRKY36, Zhongguo Nongye Kexue (Chinese Agricultural Sciences), 48(5): 851-860